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ABSTRACT 

We examine the metallicity and age of a large set of SDSS/DR6 galaxies that may be Blue Compact 
Dwarf (BCD) galaxies during quiescence (QBCDs). The individual spectra are first classified and then 
averaged to reduce noise. The metallicity inferred from emission lines (tracing ionized gas) exceeds 
by ^0.35 dex the metallicity inferred from absorption lines (tracing stars). Such a small difference 
is significant according to our error budget estimate. The same procedure was applied to a reference 
sample of BCDs, and in this case the two metallicities agree, being also consistent with the stellar 
metallicity in QBCDs. Chemical evolution models indicate that the gas metallicity of QBCDs is too 
high to be representative of the galaxy as a whole, but it can represent a small fraction of the galactic 
gas, self enriched by previous starbursts. The luminosity weighted stellar age of QBCDs spans the 
whole range between 1 and 10 Gyr, whereas it is always smaller than 1 Gyr for BCDs. Our stellar 
ages and metallicities rely on a single stellar population spectrum fitting procedure, which we have 
specifically developed for this work using the stellar library MILES. 

Subject headings: Galaxies: abundances - Galaxies: dwarf - Galaxies: evolution - Galaxies: ISM - 
Galaxies: starburst 



1. INTRODUCTION 

In the hierarchical picture of galaxy formation, large 
galaxies ar ise through the assembly of smaller aggre - 
gates (e.g.. IWhite fc Frenkl[l99lt iDiemand et alJl2007| L 
and metal-deficient dwarf galaxies are possibly the clos- 
est examples we can find of the elementary primordial 
units from which galaxies assembled. In the downsiz- 
ing paradigm, large galaxies form metals early on (e.g., 
ICowie et al.l[l9 96; Kau ffmann et al.ll2003[ ). and only low 
metallicity dwarfs may still keep a fossil record of the 
pristine Inter-Stellar Medium (ISM). For one reason or 
another, Blue Compact Dwarf (BCD) galaxies seem to 
probe early phases of the Universe. They combine 
the two required ingredients, i.e., they are dwarfs hav- 
ing record-braking low m etallicities fe.g. lKunth fe Ostlinl 
l2000t llzotov et al.ll2005T ). BCD galaxies have been used 
to constrain, e.g., the properties of the first (Pop III) 
stars that polluted the primordial ISM at the time of 
galaxy formation (e.g.. iBromm fc Larson! [2004 ) . or the 
primordial He a bundance inherited from big-bang nucle- 
osynthesis (e.g.. llzotov et al.ll2007| ). 
BCD galaxies have been ext ensively studied 

durin g the last 35 y ears (e.g. ISargent fc Searlel 

19701: iSearle fc Sartrenti Il972l iSilk et all 11987 
Davies fc Phillippd Il988t iPapaderos et alT Il996e 
levichl 119971 



Consequently, identifying precursors and descendants of 
BCDs would be a major breakthrough in deciphering the 
nature and the functioning of these special galaxies. It 
will facilitate using them as reliable tools in cosmological 
studies. 

Evolutionary connections between different dwarf 
galaxies and BCDs have been bot h proposed and 
quest io ned in th e litera tu re (e.g., ISearle fc Sargent] 
1971 ISilk et all fl987l [Davies fc Phillippsl fl98T 



Papadcro s et alj Il996at [Tellcs & T erlevichl 1 1997 
Gil de Paz fe Madorel 120051 : iDellenbusch et all 12007 
In an attempt to complement these 



Amorin et alJ 



Telles fc TerlevicI 
Gil de Paz et al.l 



2003; 

Amorin et"aH l2007f l However 



7l ICairos et al. 
iGil de Paz fc Madorel 



and evolve remains unknown, 
that characterizes BCDs lasts 



the way they grow up 
The intense starburst 
only a few Myr (e.g., 
iMas-Hesse fc Kunthl [l999t iThuanl fl99lh . BCDs seem 
to be undergoing a transient but we ignore how they 
reach such state, and what happens to them afterwards. 
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efforts, we carried out a search for galaxies that may 
be BCD during the periods where the major starburst 
is gone, i.e., quiescent BCD s or, for short, QBCDs 
([Sanchez Almeida et al.l 120081 . hereafter Paper I). We 
addressed the issue from a new perspective. Most 
BCDs show a red low surface brightness component 
which should exist before the present starburst sets 
in and should remain once the starburst fade s away 
([Loose fc Thuanl 119861 : IPapaderos et alj Il996bl ). By 
carefully removing the starburst, this underlying compo- 
nent has been studied and characterized in the literature 
(e.g- jNoeske et al.ll2003l: ICaon et aHl2005l; lAmorin et all 
[2007L I2009h . We searched the SDSS/DR6 database 
for isolated galaxies with the luminosity, color, surface 
brightness, and concentration characteristic of the low 
surface brightness component underlying the BCDs 
(Paper I). Assuming that the underlying low surface 
brightness galaxy remains unaltered after each starburst 
exhaustion, the targets thus selected could be QBCDs. 
The search yielded some 21500 QBCD candidates, with 
properties pointing out that they may be indeed pre or 
post BCDs. In particular, they have the same luminosity 
bC [4yiction as the BCDs, although they are thirty times 
more numerous. The results suggested an evolutive 
sequence where BCDs undergo many short starburst 



2 



Sanchez Almeida et al. 



phases during their life times, as proposed long ago by 
ISearle fc Sargentl (|1972| ). In between bursts, the galaxies 
show up as QBCDs in a low activity state of various 
degrees lasting thirty times longer than the bursts. 
Statistically, QBCDs should undergo a BCD phase 
every 300 Myr and lasting some 10 Myr. This sequence 
of BCD and QBCD phases can be maintained during 
several Hubble times, and the most active QBCDs 
are indeed BCDs. Paper I carries out the differential 
comparison with BCDs by selecting the sample of 
BCDs also from SDSS/DR6, and employing the same 
procedures used to retrieve the QBCDs. 

In spite of all these agreeable features, the evolutive 
link between BCDs and QBCDs presents an important 
difficulty posed in Paper I. The QBCD oxygen abun- 
dance was estimated to be 0.35 dex systematically larger 
than the oxygen abundance of the BCDs. This makes 
the role of QBCDs as BCD precursors questionable since 
starbursts increase metallicity, and the putative precur- 
sors (QBCDs) should have lower metallicity than their 
descendants (BCDs). Paper I offered a few alternatives 
to clear out the difficulty, most of which were related with 
the infall of metal poor gas before the starburst sets in. 
In addition, we speculated that the metallicity assigned 
to the QBCDs may be biased, with the true QBCD 
metallicities much lower than the observed ones, and 
close to the observed BCD metallicities. We derive the 
oxygen abundance from emission lines produced in H n 
regions, which trace the ISM in those places now going 
through a star-formation episode. In the case of QBCDs, 
the star formation rate is quite small (< 0.1 M© y _1 even 
for the brightest ones; §[8]), therefore the volume of galac- 
tic gas sampled by the measurement is very small too. 
The question arises as to whether the abundance of this 
gas is representative of the total galactic gas. If it is not, 
then it could explain a false overabundance of oxygen 
in our QBCD candidates. The sampled gas may not be 
properly mixed up with the galactic ISM and, therefore, 
be self metal enriched by successive starbursts. 

The possibility that the metallicity deduced from emis- 
sion lines may be contaminated by recent starbursts 
has been previously m e ntioned in the literature (e.g., 
Kunth & Sargentl 1 19861 iThuan et all 120041 : iDalcantonl 
200l . The mixing of the ISM is a slow process, which 
leaves behind a patchy mediu m fe.g.. lTenorio-Taglell996t 
Ide Avillez fc Mac Lowir2002ft . 

The present work was originally meant at testing the 
main conjecture in Paper I, namely, that the emission 
line derived metallicity overestimates the true average 
metallicity of the QBCD gas. If so, it should be signifi- 
cantly larger than the metallicity of other galactic com- 
ponents, in particular, the metallicity of the stars. This 
seems to be the case (§[6]), but in the way of working it 
out, several other properties of QBCDs (and BCDs) have 
emerged. These results are described here in fairly broad 
terms, keeping in mind their potential interest outside 
the specific original motivation of the work. In particu- 
lar, QBCD galaxies are quite common (one out of each 
three local dwarfs; Paper I), so that their properties may 
also be representative of the whole class of dwarf galax- 
ies. 

The paper is organized as follows: § [2] summarizes the 
main observational properties of the Sloan Digital Sky 
Survey/Data Release 6 (SDSS/DR6) spectra used in our 



analysis. § [3] explains the classification of spectra be- 
fore averaging them out to improve the signal-to-noise 
ratio. By fitting the observed QBCD spectra with syn- 
thetic spectra, we assign ages and metallicities to the 
stellar component of the galaxies (§H]). Gas metallicities 
are estimated in § [5j with their uncertainties critically 
assesses in App. EJ § H] puts forward the excess of gas 
metallicity with respect to the metallicity of the stellar 
component. Ages and stellar content are analyzed in § El 
Chemical evolution model galaxies able to account for 
the differences between stellar and nebular metallicities 
are discussed in § where we also have to estimate Star 
Formation Rates (SFRs). Finally, the main results and 
their implications are discussed in § O 

2. DATA SET: SDSS SPECTRA 

We aim at assigning metallicities to the stellar com- 
ponent of the QBCD candidates selected in Paper I. 
The original galaxies were chosen within the SDSS spec- 
troscopic catalog, to have redshifts from which we de- 
rive absolute magnitudes (mean redshift 0.030, with 
a standard deviation of 0.014). The present anal- 
ysis of stellar metallicities is based on these SDSS 
spectra. For the sake of comprehensiveness, their 
main characteristics are summarized here. A more 
detail le d account can be found in IStoughton et al.l 
(2002). lAdelman-McCarthv et all (|2008h . and also in the 
SDSS website (http://www.sdss.org/dr6l. 

The SDSS spectrograph has two independent arms, 
with a dichroic separating the blue beam and the red 
beam at 6150 A. It simultaneously renders a spectral 
range from 3800 A to 9200 A, with a spectral resolu- 
tion between 1800 and 2200. The sampling is linear 
in logarithmic wavelength, with a mean dispersion of 
l.lApix -1 in the blue and 1.8Apix _1 in the red. Re- 
peated 15 min exposure spectra are integrated to yield 
a S/N per pixel > 4 when the apparent magnitude in 
the g bandpass is 20.2. The spectrograph is fed by fibers 
which subtend about 3"on the sky. Most galaxies are 
larger than this size, therefore, the spectra sample only 
their central parts (e.g., 89% of the QBCD galaxies have 
an effective radius larger than half the fiber diameter). 
We retrieve the 21493 QBCD spectra and the 1609 BCD 
spectra in FITS format from the SDSS Data Archive 
Server. All spectra were re-sampled to a common rest- 
frame wavel ength scale that matches the spectral library 
used in § l4l (jSanchez-Blazquez et afll2006t [Cenarro et al.1 
12001 . We use linear interpolation to oversample the 
original spectra with a constant dispersion of 0.9 A pix . 
The spectra were normalized to the flux in the g color fil- 
ter (effective wavelength ~ 4825 A), a normalization fac- 
tor that we compute from each spectra using the trans- 
mission curve downloaded from the SDSS website. 

3. CLASSIFICATION OF GALAXY SPECTRA 

As we will discuss later on f§ l4.2p . the S/N of the indi- 
vidual SDSS spectra is insufficient to estimate the metal- 
licity of the stellar component. We improve the S/N ratio 
to acceptable levels by averaging similar sp ectra (a tech- 
nique often refereed to as stacking; see, e.g. lEllison et al.l 
2000). Before averaging, the spectra have been classi- 
fied in alike sets using a cluster analysis algorith m. We 
employ the simple k-means clustering (see, e.g., lEverittj 
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1 !)!)•">. Chapter 5) . A number k of template spectra are 
selected at random from the full set. Each template spec- 
trum is assumed to be a cluster center, and each spec- 
trum of the data set is assigned to the closest cluster 
center (closest in a least squares sense). Once all spectra 
in the dataset have been classified, the cluster center is 
re-computed as the average of all spectra in the cluster. 
This procedure is iterated with the new cluster centers, 
and it quits when no spectrum is re-classified in two con- 
secutive steps. The algorithm is simple and fast, but it 
yields different clusters with each random initialization - 
the final cluster centers keep some memory of the original 
randomly chosen cluster centers. This drawback does not 
interfere with our purpose of selecting sub-sets of simi- 
lar spectra suitable for averaging because, independently 
of the initialization, the clusters always contain similar 
spectra. The algorithm forces all spectra in a class to be 
similar to the cluster center, and therefore, similar among 
them. The number of clusters k is arbitrarily chosen but, 
in practice, the results are insensitive to such selection 
since only a few clusters possess a significant number of 
members, so that the rest are discarded. Figure Q] shows 
the number of elements in each class of QBCD spectra 
resulting from applying the procedure. The classes have 
been sorted and labelled according to the number of el- 
ements, with Class the most numerous, Class I the 
second most numerous, and so on (percentages are given 
in Table [TJ). Figure [5] shows the average spectrum cor- 
responding to the first nine most numerous classes. The 
classification was carried out using four spectral band- 
passes containing emission lines (the bandpasses are in- 
dicated as dotted lines in Fig [2l Class 0, and also in 
Fig. [3]). The use of these particular bandpasses empha- 
sizes the contribution of the emission lines for classifica- 
tion which, otherwise, would be completely overridden 
by the continuum. We select bandpasses throughout the 
full spectral range to assure that the global trend of the 
continuum is considered when classifying. There seems 
to be continuous variation of properties which, in the 
end, give rise to a large variety of shapes . From spectra 
without significant emission lines (Class 3), to spectra 
with red continuum and emission lines (Class 7), to blue 
continua with moderate emission lines (Class 4). The 
most numerous Class has blue continuum and presents 
emission lines. 

The scatter among the spectra belonging to a class 
depends on wavelength, and it is largest in the intense 
emission lines. As it is illustrated in Fig. [3] with Class 
spectrum, it is of the order of 10% in the spectral ranges 
with absorption lines, and it can be of the order of 50% in 
the regions having strong emission lines (see the dashed 
line in Fig. [31 which corresponds to the rms fluctuations 
among all the spectra in the class divided by the mean 
spectrum). 

The same classification procedure was also applied to 
the control set of BCDs. Representative spectra of the 
most numerous classes are shown in Fig. [4j also sorted 
and labelled according to the number of galaxies in the 
class. The most conspicuous differences with respect to 
QBCD spectra are the strength of the emission lines, 
and the (barely visible but always present) blue continua. 
The fraction of galaxies in each class is listed in Table [TJ 




10 20 30 10 20 30 

Class § Class # 



Fig. 1. — (a) Histogram with the number of QBCD galaxies 
corresponding to each class, (b) Normalized cumulative histogram, 
i.e., fraction of QBCD galaxies from Class to the each class #. 
It is given in percent. Note that the first ten classes include 90 % 
of the QBCDs. 

3.1. Green valley QBCDs 

One of the findings of the SDSS is the existence of 
a gala xy color sequence with well defined bi-modality 
(e.g., iBalogh et al1l200l . Paper I shows how QBCDs 
occupy all the color sequence between the blue and the 
red clumps. It turns out that the classification has been 
able to separate galaxies in the red sequence, galaxies 
in the blue sequence, as well as those in between (of- 
ten r efereed to as green valley galaxies, e.g.. ISalim et al.l 
2007). Figure [5] shows color vs color scatter plots for 
the QBCDs, the BCDs, as well as the most usual QBCD 
classes separately. It also includes the somewhat arbi- 
trary boundary between the red and the blue sequences 
worked out in Paper I (the dashed line). The classi- 
fication does a fair job in splitting the galaxies in col- 
ors. Among the most numerous classes, it turns out that 
only Class 3 belong to the red sequence. (Class 3 has no 
emission lines; see Fig. O) In addition, we noticed that 
Class 1 seems to include all the green valley galaxies, 
i.e., the transition galaxies, central to understand how 
and why gal axies move back and forth in the color se - 
quence (e.g.. ISpringel et al.ll2005HCattaneo et al.ll2006ft . 
The goodness of this green valley galaxy selection method 
can be appreciated by comparing the top left plot with 
the bottom right plot in Fig. El The two of them include 
the same QBCD galaxies except for Class 1. A clear gap 
splits up the red and the blue clumps. 

Since different classes have different colors, and the 
QBCD present a clear color-(nebular)metallicity rela- 
tionship (Paper I), different classes have different metal- 
licities too. Scatter plots of metallicity vs absolute mag- 
nitude are s hown in Fig. El Here an d throughout we use 
the recipe in iPettini fc Pagell (|2004f ) to compute nebular 
metallicities from the N2 strong-line rati o (see also §E1). 
It upg rades of the classic al calibration bylDenicolo et alJ 
(2002) used in Paper I ()Shi et al.ll2005l and references 
therein). Late types (e.g., Class 0) are metal poor as 
compared to the transition objects included in Class 1. 
These green valley galaxies have solar metallicity. The 
few galaxies in the red clump with emission lines (Class 7; 
see Fig. [2]) seems to have slightly super-solar metallic- 
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TABLE 1 

Properties of QBCD classes and BCD classes. 



Galaxy 


Class 


Fraction a 


Stellar 


Stellar Age c 


Nebular 


Ha EW° 


[Nil] EW f 


Type 




[%] 


Metallicity b 


[Gyr] 


Metallicity d 


[A] 


[A] 


QBCD 





36.1 


-0.44±0.03 


1.7±0.2 


-0.12±0.18 


22.3 


5.2 




1 


14.3 


-0.39±0.02 


5.2±1.4 


0.05±0.18 


5.3 


2.5 




2 


9.5 


-0.32±0.12 


1.1±0.6 


-0.25±0.18 


50.4 


6.9 




3 


8.5 


-0.33±0.02 


11.1±1.7 










1 


5.3 


-0.36±0.08 


1.0±0.1 


-0.26±0.18 


14.5 


1.9 




5 


3.5 


-0.68±0.28 


1.1±3.0 


-0.40±0.18 


99.9 


7.5 




6 


3.4 


-0.44±0.11 


1.1±0.2 


-0.25±0.18 


31.6 


4.4 




7 


3.1 


-0.42±0.04 


4.1±1.2 


0.05±0.18 


8.4 


3.9 




8 


2.5 


-0.37±0.07 


2.7±1.7 


-0.01i0.18 


27.6 


10.0 




9 


2.0 


-0.30±0.02 


17.8±0.4 








BCD 





46.9 


-0.33±0.08 


0.9±0.1 


-0.38±0.18 


81.2 


6.7 




1 


12.7 


-0.34±0.22 


1.0±0.8 


-0.41±0.18 


160.2 


11.7 




2 


7.9 


-0.37±0.40 


0.9±2.5 


-0.49±0.18 


241.6 


12.4 




3 


7.4 


-0.39±0.17 


1.1±0.8 


-0.37±0.18 


125.9 


10.7 




1 


6.5 


-0.41±0.21 


0.9±0.6 


-0.51±0.18 


142.0 


6.5 




5 


3.8 


-0.34±0.15 


1.0±0.4 


-0.35±0.18 


93.7 


8.8 




6 


3.5 


-0.34±0.54 


1.1±4.6 


-0.48±0.18 


333.0 


17.1 



Note. — It includes those classes containing 90% of the galaxies. 
a Percentage of galaxies represented by the class. b In logarithm scale, referred to the solar metallicity. 
Errors from Monte Carlo analysis in §|4] c Errors from the Monte Carlo analysis in §[4] d In logarithm 
scale, referred to the solar metallicity. Its error has been taken from Pettini & Pagcl (2004). c Equivalent 
width of Ha. No data implies line in absorption. f Equivalent width of [Nil] A6583. No data implies line 
in absorption. 



ity (Fig. E]). 

The k-mean clustering classification provides an auto- 
matic method to identify spectra of green valley galaxies. 
It works for QBCDs, however, there is no clear reason 
why it should be restricted to them. It may be valid 
for any type of galaxy, even outside the particular set of 
dwarfs we are dealing with. We elaborate on this possi- 
bility in § [9] 

4. DETERMINATION OF STELLAR 
METALLICITIES AND AGES 

The stellar content of a galaxy can be studied through 
modeling and interpretation of the absorption features 
in the integrated spectrum. The analysis of line-strength 
indic e s, mainly those of the Lick s ystem (|Worthev et al.l 
1 1994 IWorthev fe Ottavianil Il997f ). has been the most 
common approach for studying the stellar metallicities 
and ages. Most studies have been focused on early- type 
galaxies, wi th stellar populati ons typically older than 
1 Gyr (e.g., iTrager et al.lll998l and references therein). 
The use of age-sensitive Balmer line indices, e.g., H/3, and 
metallicity-sensitive indices, e.g., Mg6, allows us to lift in 
part the age-metallicity degeneracy affec ting the stellar 
popu l ations of early-type galaxi es (e.g., IWorthev et alJ 
1 1994 IVazdekis fc Arimotol Tl999f ) . However, these line 
strengths, and particularly the Balmer indices, are not 
optimal for analyzing our galaxy spectra, since they are 
filled-in with nebular emission (Fig. [2|). Here we take a 
different approach, and ages and metallicities are esti- 
mated by direct comparison of the observed spectra with 
model spectra from stellar population syntheses. The full 
spectral range is used simultaneously. Th is alternative 
strategy has a respectable tradition (see iKoleva et al.l 



120081 and references therein), and it allows us to eas- 
ily overcome the problem of emission lines by masking 
them out. Emission lines represent only a small frac- 
tion of the spectral range, and the rest of the spectrum 
can be used to extract the required information (sec 
below). We us e an updated version of the models by 
iVazdekia (|1999l) , which provide spectral energy distribu- 
tions of single-age, single-metallicity stellar populations 
SSPs) on the basis of the stel l ar spectral library M ILES 
Sanchez-Blazquez et al.ll2006HCenarro et al.ll2007f l. The 
MILES SSP spectra combine, according to a Salpeter 
initial mass function distribution, a suite of stellar spec- 
tra from 0.09 M Q to 100 M Q . They have a resolution of 
~ 2.3 A, a spectral range from 3540 A to 7410 A, and a 
dispersion of 0.9 A pix . MILES SSP extends the range 
of ages ofVazdekis ( 1999|, and now it covers from 0.1 Gyr 
to 17.8 Gyr. MILES SSP spectra span a range of metal- 
licities 2 between —1.7 and +0.2. The grid includes 276 
SSP spectra, with 46 samples equispaced in logarithmic 
time, and 6 steps in metallicity. The range of metallic- 
ities and ages fits in well the values to be expected for 
QBCD (see §[J). As far as the wavelength sampling and 
wavelength coverage are concerned, MILES SSP spec- 
tral resolution is comparable to SDSS (although better), 
but it misses the reddest 1800 A of the SDSS spectral 
range. The uncovered 20% of the SDSS spectral range 
is in the near IR, where the number of spectral lines de- 
creases significantly. Keeping in mind all these variables, 
MILES SSP meets very well our needs. 
The fits are carried out by direct comparison of each 

2 In the usual logarithmic scale refereed to the solar metallicity 
Zq, i.e., log(Z s /ZQ) with Z s the fraction of mass in metals. 
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Fig. 2. — Average spectra of the first nine most abundant QBCD classes. Together with the class identifier, the insets of the figures give 
the number of galaxies in the class. Wavelengths are in /im, and the range of ordinates is different for the different plots. All spectra are 
normalized to the flux in the g filter. The dotted line shown together with Class indicates the band-passes used to classify spectra (i.e., 
those wavelengths were it is zero were disregarded for classification). 



average profile representative of a class with all spectra 
in the MILES SSP library, smeared to three spectral res- 
olutions (the original one, the original one plus 2.5 A, 
and the original one plus 3.5 A). Considering various 
broadenings is required to account for stellar motions, 
as well as for the difference of spectral resolution be- 
tween MILES SSP and SDSS. The observed spectrum is 
compared with each synthetic spectrum, and closest one 
in a lest-squares sense is chosen as best fit. The com- 
parison was carried out with a few constraints which try 
to minimize potential biases. (1) A 100 A running-box 
mean of the original spectra was removed from observed 
and synthetic spectra. By removing the continua, the 
results of the fits are not very sensitive to the extinction, 
a miscalibration of the spectra, the uncorrected differ- 
ential refraction (|Izotov et all [20061 . and so on, which 
affect the continua but not so much the relative intensity 
of adjaccent spectral lines. Moreover, it guarantees that 
ages and metallicities are inferred from spectral lines, 
with negligible contribution from the global shape of the 
continuum. (2) We assume the observed spectra Oi to 
be a linear combination of a starburst spectrum N plus 
an stellar spectrum Si, 

O l =aN l +(3S ll (1) 

with the underscript i representing the i — th wavelength, 



and a and (5 being two scaling constants. The starburst 
spectrum has strong emission lines and little continuum, 
therefore, one could simply neglect the core of the emis- 
sion lines when carrying out the fits (i.e., one could mask 
out the emission lines and set Ni = in equation Q]). 
Here we go a step further so that the (small) contamina- 
tion by Ni outside emission lines is estimated and sub- 
tracted out. The decontamination procedure works as 
follows. At the emission line cores of Oi, the spectrum is 
dominated by JVj so that, 

O l (1-Wi) = (l-Wi)(aNi+0Si) si (l- Wi )aNi, (2) 

with Wi a properly chosen weight which is zero in the 
emission cores and one elsewhere, i.e., 

w _ f emission lines, /g\ 
1 \ 1 elsewhere. 

Using the different classes of BCD spectra as proxies for 
Ni, we choose for each Oi the N BCD spectrum that 
minimizes the appropriate merit function, 

X^^iiOi-aNifO-Wi) 2 ], (4) 

i 

with 

a = [ J2 Oi N (1 - m ( ) 2 ] /E^K 1 - w 3 f] , (5) 

i 3 
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Fig. 3. — Expanded version of the spectrum of Class in Fig. [2] In addition to the spectrum itself (the solid line), the plot includes the 
rms fluctuations among all the spectra in the class divided by the mean spectrum (the dashed line). The bandpasses used for classification 
are shown as dotted line. The plot also labels several typical spectral features. W avelengths are given in A. The large fluctuations between 
6200A and 6300A are produced by telluric lines (e.g., § 4.2 in[Gravcs et al. 2009). Their presence does not affect our fits (see the residuals 
at these wavelengths in Figs. [7]and[8jl. 



the latter being just a least squares estimate of the scal- 
ing factor that best fit the emission lines of Oi once TV; is 
given. Note that the weight (1 — to*) in equation (01 as- 
sures that only the emission lines contribute to \ 2 ■ The 
Ni and a thus derived allows us to compute the observed 
spectrum corrected for emission, O*, 

Of = Oi - aNi = 0Si. (6) 

The best fitting MILES SSP spectrum is obtained by 
repeating the same procedure with O* but masking out 
the emission lines, i.e., defining the merit function for 
each MILES SSP spectrum S t {t, Z s ) as, 

X 2 (t, Z.) = [{OI - PS % {t, Z s )) 2 w 2 } , (7) 

i 

where 

P= [E O* S i {t,Z s )w1]/[Y,S j {t,Z s f w^]. (8) 

i J 

The expressions explicitly include the dependence of the 
synthetic spectrum on the age of the starburst, t, and 
the stellar metallicity, Z s , i.e., Si(t,Z s ). The weight if, 
in equation JJJ) cancels out the contribution of the emis- 
sion lines, rendering the correction ([5]) of secondary im- 
portance. The weights iWj are assigned so as to cover 



the emission line cores observed in BCD spectra. Exam- 
ples of these weights are the (thin) solid lines in Figs. [7] 
and [5] The positions of the minima of these broken lines 
mark the wavelengths discarded from the fits. The fit- 
ting procedure described above was also applied to BCD 
spectra. In this case we cannot correct for the starburst 
since BCDs are used as template starburst spectra. We 
just mask out the emission lines and force a = in equa- 
tion ©. 

As judged from visual inspection, the best fitting 
model spectrum reproduces very well the observed spec- 
tra; see, e.g., Fig. [7] Error bars cannot be assigned us- 
ing traditional methods based on the Hessian matrix of 
X 2 , since we do not have a continuous function x 2 (i, Z s ) 
to compute the partial derivatives with respect to age 
and metallicity. We resort to Monte Carlo simulations 
to assign confidence intervals. Gaussian noise with the 
standard deviation of the residuals is added to the best 
fitting model spectrum. This mock observation is ana- 
lyzed as the real observation to get an age and a metal- 
licity which, in general, differ from those of the best fit- 
ting model spectrum. The procedure is repeated 1000 
times, which provides a range of ages and metallicities 
consistent with the best fitting model spectrum and the 
residual of the fit. Confidence intervals thus assigned re- 
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Fig. 4. — Same as Fig. E] but for BCD galaxies. 



veal a serious problem of degeneracy in the metallicity 
estimate. The standard deviation of the most common 
QBCD Class turns out to be 0.25 dex, which allows for 
any metallicity between —0.20 and —0.70. Such degener- 
acy in metallicity is the young stellar population equiva- 
lent of the well known age-metallicit y degeneracy ap pear- 
ing in early-type galaxy dating (e.g.jWorthcy 1994). We 
managed to sort out the degeneracy problem by over- 
weighting the contribution of those spectral bandpasses 
that are known to be particularly sensitive to metallic- 
ity. Specifically, the merit function in the definition Q 
is now, 

{0 emission lines, 
W metallicity sensitive bandpasses, (9) 
1 elsewhere, 

with W > 1 for overweighting. Th e band passes were se - 
lected from the Lick index system (jWorthev et al.|[l99^ 1. 
which was specifically designed for estimating ages and 
metallicities in the integrated light of stellar populations. 
In order to find out which are the Lick indexes most 
sensitive to metallicity in our domain of ages, we com- 
puted the variation of the indexes with metallicity at 
a given age. Some results for the MILES SSP library 
are shown in Fig. [5) Among the 21 indexes defined by 
iWorthev et al.l (|1994l ). we select the three indexes in the 
top row because they show the largest variation with 
metallicity. Figure [9j bottom row, also includes three 



other indexes commonly used in metallicity studies of 
early-type galaxies (e.g., t hey are combined to form the 
so-ca lled [MgFe] index; see iGonz alcz 1993; lThomas et all 
l200"l . The range of variation is clearly inferior to the 
variation of the indexes that we select. The bandpasses of 
the three selected indexes, Fe4383, Fe4531, and Fe4668, 
are indicated in Figs[7K and UK as the wavelengths where 
the thin solid line representing the weight w; goes out of 
scale. We tried with various overweights (W = 10, 20, 50 
and 100), to finally choose W = 50 since the trial fits 
indicate that the inferred metallicity does not depend 
on the actual weight when the weights are large enough. 
The use of these overweights improves the metallicity es- 
timate to a large extent. We repeated the Monte Carlo 
simulation described above, and the random errors of 
QBCD Class decrease by almost an order of magni- 
tude with respect to the case where W was set to one. 
Table [T] includes the standard deviation for the ages and 
metallicities of all major QBCD and BCD classes. They 
will be employed as 1 a errors in the discusions along 
the paper. The small value of these random errors has 
been indepen dently corroborated by a bootstrap error 
estimate fe.g.. lMoore et al.ll2003t ). A caveat is in order, 
though. These small errors only indicate that the best 
fitting MILES SSP spectrum is well defined, i.e., among 
the MILES SSP set, a few spectra reproduce the obser- 
vation clearly better than the rest. Our procedure does 
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Fig. 5. — Color vs color scatter plots for the QBCDs, the BCDs, and the most usual QBCD classes separately. The insets specify each 
galaxy group. The dashed line is the same in all plots, and it was worked out in Paper I to separate the blue and the red sequences. Note 
that most QBCDs belong to the blue sequence. Almost all galaxies in the red clump correspond to Class 3. Moreover, Class 1 seems to 
gather most galaxies in the so-called green valley between the red and the blue sequences (c.f. the plots in top left and the bottom right 
corners, which are identical except that Class 1 has been removed from the latter). 
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Fig. 6. — Scatter plot of oxygen abundance vs g—r color for some 
representative QBCD classes. The horizontal solid l i ne co rresponds 
to the solar metallicity as given by Grcvcsse ct al. (2007). Class 3 
is not included because it lacks of the emission lines needed to 
compute nebular abundances. 



not account for systematic errors, which may dominate 
the error budget (e.g., is SSP a good description of our 
galaxies?). The magnitude of the systematic errors is 



unknown, and ignored in our discussions. 

4.1. Self- consistency of the continua 

A running mean average was subtracted from both the 
observed and the model spectra to minimize the influence 
of miscalibrations (§S]). Then our fits are virtually blind 
to the galaxy continua. The question arises as whether 
the ages and metallicities thus derived are or not consis- 
tent with the observed galaxy continua. 

In order to assign a continuum to the model spectra, 
one has to bring out the continuum information removed 
when subtracting the running mean average. We do it 
by parameterizing the relationship between the actual 
observed spectrum, Oi, and model we fit, mj, including 
the biases that the subtraction of a continuum removes, 



1Q - (Ai -A )/2. 5 + 



(10) 



where Ai co rresponds to extinct ion by dust 3 defined as 
usual (e.g., ICardelli et al.1 11989[ ), and k, accounts for 
other possible differences not included in the model. (Re- 
call that the underscript i parameterizes the variation 

3 Both internal, and due to our Galaxy, since the two of them 
add up when dealing with low redshift targets. 
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Fig. 7. — Observed average spectrum of QBCD Class (the thick solid line) and best fitting MILES SSP synthetic spectrum (the dotted 
line). The dashed line corresponds to the residuals vertically shifted by an arbitrary amount. Panels (a) and (b) show the full spectral 
range, whereas (c), (d), (e) and (f) zoom into details to appreciate the goodness of the fit. The weights of the fits are represented in its 
own scale as a thin solid line, with the minima corresponding to no contribution, i.e., to weight equals zero. The weight goes out of scale in 
(a) indicating the wavelengths of the three Lick metallic indexes overweighted during fitting (Fe4383, Fe4531, and Fe4668). The continuum 
has been subtracted from both the observed, and the synthetic spectra. 



with wavelength.) Equation (|10[) also assumes that the 
model spectrum includes a gray extinction given by Aq 
(incorporated into the global scaling factor (3 that we 
use for fitting; see equation [7]). As we will show, the ex- 
pression (TlTJ|) is fully consistent with equation ([T]) if one 
removes a running mean average of the observed spec- 
trum, (o,-). The spectrum we fit (Oi in equation [T]) is 

O t = Oi - (oj> ~ (aN t + (3Si) lQ-(^-^)/2. 5) (n) 

with the model galaxy spectrum given by 

ceNi + pSi = mi-(mi). (12) 

We have employed equation (|10[) assuming that Ai and 
Ki do not vary within the kernel that defines the running 
mean (i.e., (Ai) — Ai, and = Ki). Neglecting in 
equation (fTTj) terms of the order of, 

(aNi+0Si)(Ai-Ao), (13) 
one ends up with, 

Oi~aNi + pSi, (14) 

which is the approximation used for fitting (equa- 
tion Q]). Within this approximation, one can re- write 
equation (fTU|) as, 

o l ~aN l +(3S l + (o l ), (15) 



where 

(<*) = [1 - {Ai - A )/(2.5 loge)] (to,) + K h (16) 

is the term to be added to the best fitting synthetic spec- 
trum, aNi + (3 Si, to recover the observed spectrum with 
continuum, o\. Equation (|16|) allows us to estimate both 
Ai and Ki and, consequently, to complete the best fitting 
model with its continuum. For lack of better assumption, 
we regard Ki as independent of wavelength. In addition, 
the wavelength dependence of Ai is assumed to be known 
except for a scaling factor, parameterized as the extinc- 
tion in the Johnsson's V band Ay. Th e ratio Aj/Ay 
is assu med to follow the milk y- way law bvlCardelli et al.l 
(|1989fh modified according to lMisselt et al.l(|1999D to rep- 
resent the large Magellanic cloud, which we use as a 
proxy for low metallicity extinction law. (The conclu- 
sions below remain even if one directly takes the milky- 
way extinction law.) Then the constants Ki and Ay can 
be retrieved from a linear least squares fit using equa- 
tions (fl5|) and (fl"6]) since Oj, m», (rrii) and Ai/Ay are all 
known, and one can regard Aq as the (wavelength) av- 
erage extinction. The comparison between observations 
and model spectra including continuum for the nine first 
QBCD classes is shown in Fig. [TUJ (The emission lines 
have been artificially taken out to better appreciate dif- 
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Fig. 8. — Same as Fig. [7] except that the spectrum corresponds to Class 3 QBCD, i.e., the class corresponding to the red sequence 
galaxies without significant emission lines. 
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Fig. 9. — Variation with metallicity of various Lick indexes in 
MILES SSP spectra. Each curve of each plot corresponds to a 
constant age. We only show young populations, with ages between 
0.5 Gyr (the curves of smallest equivalent widths) and 2 Gyr (the 
curves of largest equivalent widths). The wavelength range cor- 
responding to the indexes in the top row are overweighted in our 
fits to break down the metallicity degeneracy. Other commonly 
used indexes are discarded because they present less dependence 
on metallicity for this range of ages; see the bottom row. The index 
name can be found in the ordinate axis labels. The corresponding 
bandpasses are defined in Worthey et al. (1994). Metallicities are 
given in a logarithm scale referred to the solar metallicity. 



ferences between observed and model continua.) Ay is 
forced to be non-negative, so that if a (small) negative 
number is found in an unconstrained fit, it is automati- 
cally set to zero. The agreement is good, in particular, 
for the most numerous classes. Keep in mind that the 
fitting procedure disregards continua, yet, the observed 
and model continua match quite well. The agreement is 
found for low extinctions, of only a few tenths of magni- 
tudes, Ay = 0.18 ± 0.27. Moreover, the most populated 
classes are in the low extinction range of such interval, 
e.g., A v = 0.04 for Class (see the labels in Fig. [T0|) . 

4.2. Why we do not use spectra of individual galaxies 
to estimate ages and metallicities 

Only average spectra are used in our analysis. In- 
sufficient signal-to-noise ratio refrain us from assigning 
ages and metallicities to individual galaxies. The reason 
stands out clearly from the error budget analysis in § 21 
The rms fluctuations of the residual of the fits are as 
small as 1.5% (see Figs. [7] and [8]) and, even in this case, 
the constraint they provide are quite loose. The individ- 
ual SDSS spectra have S/N > 4 , and this sole random 
error would rise the residual of any fit to an rms < 25%. 
This residual is some 15 times larger than the residuals 
of our fits, and such a large error would make our analy- 
sis completely unreliable. One spectrum is not sufficient. 
Putting the same idea in other words; if the errors of the 
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Fig. 10. — Similar to Fig. [2] except for the ordinate scale, magnified to appreciate differences between the observed continua (the solid 
lines), and the one inferred from fitting absorption lines (the dotted lines). The emission lines have been taken out to avoid overcrowding 
of the plots. The insets provide the class number of the spectrum together with the extinction coefficient Ay- 



stacked spectra are similar and independent, averaging 
at least some ~ 15 2 ~ 220 spectra is required to get the 
kind of residual represented in Figs. [7] and [5] 

5. GAS METALLICITIES 

The gas (or nebular) metallicities of the different 
classes of QBCDs and BCDs are estimated using 
strong line empirical calibration methods. Specifically, 
we primarily use the so -called N2 method as provided by 
iPettini fc Pagell (|2004f h We cannot employ the more ac- 
curate T e estimate because the [OIIIJA4363 line required 
to co mpute electron temperatures (see, e.g.. Hzotov et al.l 
120061 § 3.1) is much too faint in QBCDs. bmce us- 
ing strong line methods is always c ontroversial (e.g., 
lStasihskall200l [20081: iShi et all [20051 ) . we have studied 
some of the potential biases that may arise. The suc- 
cess of empirical calibrations resides in the agreement 
between the physical conditions of the calibration tar- 
gets, and those of the galaxies to be analyzed. As we dis - 
cuss in App.lAl the calibration by IPettini fc Pagell (|2004l ) 
holds for a fairly large range of conditions, broad enough 
as to encompass the different physical conditions to be 
expected in QBCDs and BCDs. Moreover, the various 
available strong line calibrations give consistent results 
when applied to our spectra. If the observed differences 
between the metallicities of BCDs and QBCDs (Paper I, 
but also § [1] and the forthcoming paragraphs) were an 



artifact of using strong line calibrations, different cali- 
brations should provide different biases. However, they 
coherently show the QBCDs to be more metallic that the 
BCDs. Figure [TT] prese nts estimates ba s ed on N2 and 
03N2 as calibrated by IPettini fc Pagel (|2004l ). When 
applied to our spectra, the two of them agree within 0.1 
dex; compare the squares (N2) and the asterisks (03N2) 
in Fig. [TTJ We have also tried with the S 23 index as cal- 
ibrated by iDfaz fc Perez- Monterq (|2000l ). giving results 
similar to N2 and 03N2 4 . In short, the difference of gas 
metallicity between QBCDs and BCDs does not seem to 
be a bias caused by using the N2 method. We do not 
correct for dust extinction to derive metallicities. This 
approximation can be readily justified since our main 
calibration N2 uses two spectral lines so close in wave- 
length that the correction for extinction is truly negligi- 
ble (~ 0.0006 dex for a one magnitude extinction). We 
measure the mean ext i nction for QBCDs to very small 
(§ [LI]), and IWu etafl (|2008fl show how the metallicity 
measured in a few representative BCD starbursts is not 
biased by extinction. 

The gas metallicities thus obtained are absolute - in 
the end, they are based on photo-ionization modeling 

4 Other methods, like P23 and P (IShi et al.l 120051) . cannot be 
applied because some of the required emission lines lie outside the 
spectral range of the SDSS spectra. 
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which relates the number of observed photons with the 
numb er of emitting atoms in the photoionizcd nebula 
(e.g., [Stasihska 2004]). The stellar metallicity, however, 
is relative to the solar metallicity. In order to com- 
pare gas and stars, the gas metallicity must be normal- 
ized to the solar value. This normalization is delicate 
and may bias the comparison, particularly in this mo- 
ment when a major revision of the solar metallicity scale 
has occurred (lAilende Prieto et al.ll200H I Asplundl 120051 ; 
iGrevesse et al.l 120071) . and it is not consistently imple- 
mented in modeling. For this reason, we feel compelled 
to discuss our use of the modern oxygen abundance for 
gas metallicity normalization, 

12 + log(0/H) o = 8.66 ±0.05, (17) 

despite the fact that the MILES SSP spectra used 
in our stellar metallicity estimates are based on a 
stellar librar y whose metallicities date back to pre- 
revision davs (ICenarro et al. 1120071 1 Le Borgne et all2004t 
ICavrel de Strobel et al.ll2001l ). The modification of the 
of solar metallicity had to do with improved modelling 
- NLTE effects and realistic 3D hydrodynamical model 
atmospheres h ave been incorporated into the analysis 
(IAsplundll2005h . Since the observed solar spectrum has 
not been modified, the revision simply re-labelled it with 
a different metallicity For the sake of argumentation, as- 
sume that the spectrum of one of our galaxies has solar 
metallicity. Then it has the metallicity corresponding to 
the oxygen abundance in equation (jTTJ) , rather than the 
metallicity originally assigned to it. Consequently, the 
use of MILES SSP spectra to estimate stellar metallici- 
ties is consistent with using equation (JXTJ) for the solar 
metallicity that normalizes the absolute gas metallicity 
inferred from emission lines. Table Q] lists the relative N2 
gas metallicities thus computed, together with the equiv- 
alent widths of two lines used in such estimate (Ha and 
[Nn] A6583). 

6. GAS METALLICITY VS STELLAR 
METALLICITY 

In this section we compare the stellar metallicities 
worked out in § 0] with the gas metallicity derived from 
emission lines in § O 

The scatter plot gas metallicity vs stellar metallicity 
is shown in Fig. [TlJ It includes all QBCD classes except 
Class 3, which has no emission lines. Several features are 
notable. First, the stellar metallicities of the most rep- 
resentative Classes (0 and 1) are systematically smaller 
than the nebular metallicities 5 . This result holds even 
when the (large) error bars of our metallicity estimates 
are taken into account. Figure [TT1 includes the error bars 
assigned to the metallicities of Class 0, and it is clear that 
the retrieved stellar and nebular metallicities disagree. 
The stellar metallicity error bar has been taken from the 
Monte Carlo simulations described in § |4] and is listed 
in Table [TJ As for t he emission line metal licity, we take 
0.18 dex inferred bv lPettini fc Pagell (|2004fl from the dis- 
persion of the N2 based metallicity when compared with 

5 Since these metallicities are refereed to the solar value, we can 
directly compare the global metallicity by mass provided by the 
spectral fitting procedure, with the oxygen metallicity by number 
inferred from emission lines. They correspond to the same quantity 
under the implicit assumption that the relative metal abundance 
of our targets follow the solar composit ion. This should be a good 
approximation for dwarf galaxies (e.g., Michiclscn ct al. 2008). 




Stellar Metallicity 

Fig. 11. — Emission line based metallicity vs stellar metallicity 
for the set of QBCD classes. The size of the symbol indicates the 
number of galaxies represented by the class, as specified in the 
inset. Boxes and asterisks correspond to two different estimates 
of emission line oxygen metallicity. The numbers inside the square 
symbols identify the major classes. Error bars for Class are shown 
for reference, and they are similar to those of all major classes; see 
Table [T] 

the more precise T e method. (See also App. |A]) Error 
bars are similar for other classes. They have not been in- 
cluded to avoid cluttering Fig. [Ill but they are listed in 
Table[TJ Note that contrarily to the behavior of Classes 
and 1, Classes 2 and 4 present similar stellar and neb- 
ular metallicities. This is not so much due to a change 
of stellar metallicity, but to a significant decrease of the 
nebular metallicity. The same kind of agreement at low 
metallicity occurs for all BCD classes. Figure fT2l includes 
the scatter plot of nebular vs stellar metallicity for the 
BCD classes (the solid star symbols). Stellar and nebular 
metallicities agree in this case, discarding serious system- 
atic errors biasing our conclusion. To be more precise, 
assuming that systematic errors in QBCDs and BCDs 
metallicities are similar, the QBCDs and the BCDs tend 
to have the same stellar metallicity but different nebular 
metallicities. Taking the most numerous Class to rep- 
resent them (i.e., the largest symbols with error bars in 
Fig- EH), the stellar metallicities Z s of QBCDs and BCDs 
are similar, 

log(Z s /Z Q ) QB CD log(Z s /Z ) BC D, (18) 
and also similar to the nebular metallicity of BCDs, 

[0/H] BC D^l0g(Z s /Z Q )BCD, (19) 

which differs from the nebular metallicity of QBCDs, 

[0/H] QBCD - [0/H] BC d + 0.35. (20) 

As usual, we have employed the notation where [O /H] = 
log(0/H) — log(O/H) , and Zq stands for the solar 
metallicity. Equations (fT5)) . ([HI)) , and (|2H)) combined 
yield, 

[0/H] QBC d =a log(Z s /Z ) QB cD + 0.35. (21) 

All the above identities have an uncertainty of the order 
of 0.2 dex, which is large but does not invalidate the 
trends. 

Figure[T2l the dashed line, includes the relationship be- 
tween nebular and stellar metallicities corresponding to a 
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Stellar Metallicity 

Fig. 12. — Same as Fig. [TT] but for both QBCDs (the square 
symbols) and BCDs (the star symbols). The dashed line shows the 
relationship for star-fo rmi ng galaxies inferred fr om the works of 
Trcmonti ct al. (2004) and Gallazzi ct al. (2005) (see main text). 
The size of the symbol scales with the number of galaxies in the 
class. Our estimate of Class error bars is shown for reference. 



large set of SDSS star-forming galaxies. We have inferred 
such a relationship by combini ng the medians of the 
mas s-metallicity re l ation sh ips by ITremonti et al.l (|2004l ) 
and lGallazzi et all (|2005f) . ITremonti et all (120041) derive 
metal licities from emission lines, whereas Gallazzi ct al. 
(2005) metallicities refer to the luminosity weighted 
stellar metallicity. The emission line metallicities are 
referred to the solar met allicity using the solar oxy - 
gen abundance employed bv lCharlot fc Longhetti| (120011) , 
which is the code behind ITremonti et al.l (|2004f ) esti- 
mates. Figure shows how the nebular metallicities 
are systematically larger than the stellar metallicities for 
galaxies of the same mass, and such difference is similar 
to the one we find for QBCDs (equation [H]). We cite 
this disagreement as a consis tency test for our me tallic- 
ity determinat i ons si nce both ITremonti et alj (|2004f ) and 
IGallazzi et al.l (|2005l ) derive metallicities using tools dif- 
ferent and more elaborated than the ones used here. If 
an unknown bias is causing the differences between neb- 
ular and stellar metallicities in QBCDs, it does not seem 
to be due to our specific simplifying hypotheses. 

The error bars employed so far correspond to la, or 
68% confidence level. If we use 2 a instead, we cannot 
discard the agreement between the gas and the stellar 
metallicities of QBCDs (in the c ase of N2 based nebu - 
lar metallicities, 2 a ~ 0.41 dex; iPettini fc Pagell2004l) . 
Note, however, that the error bars used for the nebu- 
lar metallicity are rather conservative. We find that the 
metallicities inferred from 03N2 and N2 agree (§ [5] and 
Fi g.nTI), and the sc a tter o f the 03N2 relationship found 
by IPettini fc Page! (|2004f ) is significantly smaller than 
that of N2 method ( 2 a ~ 0.25 dex fo r 03N2 ). Moreover, 
the scatter found by IPettini fc P agcl (2003) corresponds 
to individual extragalactic H II regions. Part of such 
scatter have to be of random nature, and it cancels when 
averaging many different regions or, as we do, many dif- 
ferent galaxies. Only the (unknown) systematic part of 
the error would be of relevance in our case. 

7. AGE OF THE STELLAR COMPONENT 








1 




1 1 1 


llicity 


-0.5 




a □ 




a) 

a 










in 


-1.0 










-1.5 


1 




QBCD □ - 
BCD * 

, , l 







1 1.0 




10.0 



Stellar Age [By] 

Fig. 13. — Stellar metallicity vs stellar age scatter plot for both 
QBCDs (the square symbols) and BCDs (the solid star symbols). 
Error bars for Class QBCDs and BCDs are included. The size 
of t he s ymbols codes the number of galaxies in the class, like in 
Fig. 1111 The range of the axes corresponds to the full range of ages 
and metallicities spanned by the MILES SSP library. Note how 
the assigned ages and metallicities occupy a well-defined narrow 
region among of the possible solutions. 

The absorption line spectrum fitting procedure in § [4] 
provides ages and metallicities for the stellar component 
of the galaxies. Figure Q1J] shows the metallicity vs age 
scatter plot corresponding to the two sets of galaxies, 
QBCDs and BCDs. We find that QBCDs are systemat- 
ically older than BCDs. QBCDs have ages in excess of 
1 Gyr whereas all BCDs have ages inferior to (but close 
to) 1 Gyr. The case of QBCD Class 3 is worthwhile 
mentioning separately (the oldest age in Fig. [T31 of the 
order of 11 Gyr). It corresponds to the QBCD galax- 
ies without emission lines (Fig. [2]), which form the red 
clump of the color sequence (Fig. [5|). These properties 
hint at Class 3 being early type galaxies, and the age we 
find is also consistent with this possibility. Such a very 
old origin of the stellar population of Class 3 is very well 
constrained according to our error analysis - see Table [TJ 

According to the conjecture we are examining in the 
paper (§[J), QBCDs undergo successive starbursts that 
transform then to BC Ds for short periods (lastin g only 
10 Myr or so; see, e.g.. iMas-Hesse fc Kunth)ll999[ ). The 
number density of BCDs and QBCDs requires the burst- 
ing phase to appear, statistically, every 0.3 Gyr. The fact 
that this timescale differs from the age we assign to the 
QBCDs is not at odds with the conjecture. The stellar 
population we observe now has been produced not just 
during the last starburst, but during several bursts. It is 
the luminosity weighted age of these populations what we 
have estimated, which has to exceed the age of the last 
starburst. The fact that the age of the stellar population 
of BCDs is shorter but not very different from the age of 
QBCDs adds on to this picture. If BCDs and QBCDs are 
basically the same galaxies, but the BCDs happen to be 
in a phase of enhanced star formation activity, then the 
underlaying stellar populations must have similar proper- 
ties. This turns out to be the case. The metallicities are 
similar (§[6] and equation [18]). and BCD stellar ages are 
shortened due to the strength of the current starburst. 

The fact that the age of the stellar population of BCDs 
is much larger than the age of a typical starburst sup- 
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ports that BCD are not forming stars for the first time. 
They have an underlaying stellar population much older 
than a star burst. 

8. METAL ENRICHMENT AND STAR FORMATION 

This section analyzes the difference between the ob- 
served nebular and stellar metallicities of QBCDs. We 
will find that the nebular metallicity of QBCDs seems to 
be much too high to be representative of the galaxy as a 
whole. Should it be representative, the QBCD galaxies 
would have been forming stars during the last few Gyr 
at an unobservedly large SFR. As we conjectured in Pa- 
per I, the metallicity of QBCDs inferred from emission 
lines probably represents a small fraction of the galactic 
gas, locally contaminated by recent starbursts. 

In principle, the difference between stellar metallicity 
and gas metallicity found in § [6] may be explained in 
terms of the chemical enrichment of the ISM during the 
time span between the formation of the stars and the 
present epoch. Assume that the chemical enrichment 
has followed a closed-box evolution. (The consequences 
of an open-box evolution will be discussed later on.) In 
this case the conservation of metals imposes the following 
constraint between the mass of stars M a , the mass of gas 
M g , the metallicity of the stars Z s , the metallicity of the 
gas Z, and the yield y, 



Z M g + Z s M s = y M s + Z M. 



(22) 



(We will use without explicit citation well known results 
from the theory o f chemical ev olution of galaxies; see, 
e.g. iTinslevi [1985 iPagell 119971 ) The left hand side of 
equation (|22p gives the amount of metals now existing 
in stars and gas, which is equal to the metals created 
by stars, y M s , plus the metals existing at the beginning 
of the starburst, Zq M, where M stands for the total 
mass in the star-forming closed-box, and Zq represents 
the initial metallicity of the gas. Equation ([22]) can be 
rewritten in a more convenient form, 



M = 



y + Z Q - Z s 



(23) 



y + Z-Zs ' 
with n the mass fraction of gas, 

,i = Mg/M. (24) 
Also from the theory of closed-box evolution, 

Z — Z = — yhifj,, (25) 

and 



A* = 1- 



(1 - R) SFRt 
M ' 



(26) 



with SFR the average Star Formation Rate during the 
past time interval t, and R the fraction of stellar mass 
that returns to the ISM rather than being locked into 
stars and stellar remants. Equations ([23]) and ([25]) com- 
bined provide the difference of metallicity between gas 
and stars, 

Z-Z s = -y[^ + l}. (27) 
L 1 — fi J 

In principle, the mass of the starburst M is unknown, 
but our ignorance can be parameterized using a scaling 
factor / between M and the mass of stars in the galaxy 
at present M*, i.e., 

M = /M„. (28) 



Although it is not a primary observable, the stellar mass 
content of a galaxy can be inferred from its observed 
lu minosity and co l or. U sing the calibration modeled 
bv lBell fc de Jonel (|2001h . the color transfo rmations be- 
tween Johnsons's colors and SDSS colors bv I Jester et al.l 
(2005), and the typical colors of the QBCD candidates 
in Paper I, one finds a relationship between the stellar 
mass in solar mass units, and the absolute magnitude in 
the SDSS g color, 



log(M*/M ) 



-0.50 g+ 0.35. 



(29) 



The parameters R and y are constrained by the stellar 
evolution models, so that R ~ 0.2 ( e.g., ITinslevi 119801 : 



IPagell 119971 : lAparicio fc Gallartl l2004h and, for oxygen, 
y ~ 4 x 10~ 3 (|Dalcantonll2007l and references therein). 
Then given the age, the SFR, and the mass (i.e., /) of a 
starburst, equations ([26]) . ([27]) and ([28]) allows us to pre- 
dict the difference of metallicity between gas and stars, 
which is the parameter measured in § [S] Consequently, 
we can use them to estimate the mass of the starburst 
and/or the SFR required to explain the observed metal 
enrichment of the gas. This is what we do next. 

One can estimate the present SFR of the QBCD galax- 
ies fro m their Ha emiss ion. We have used the prescrip- 
tion in iKennicutt] |l998), which gives the SFR as a func- 
tion of the Ha luminosity. We derive the luminosity from 
the observed Ha equivalent width, and the absolute lu- 
minosity of the galaxies in the SDSS r bandpass, approx- 
imately centered at the Ha wavelength. The transforma- 
tion between magnitudes and fluxes has been carried out 
keeping in mind that the SDSS color system is an AB 
system ([Smith et al.l l2002). which renders, 



SFR: 



7 r i( J 

100 A 



■0.4(r+19.0) 



M Q yr- 



(30) 



with W}i a the equivalent width in A, r the integrated 
absolute r magnitude, and 7 the fraction of galactic 
light contributing to the the starburst. 6 Figure [T4l shows 
the SFRs inferred from equation (|30p for the individual 
QBCD galaxies in Class (the square symbols). We 
have used 7 = 0.5, as a reasonable upper limit for the 
extent of the starburst, but using 0.1 or 1 do not change 
any of the conclusions discused below. Typically, the 
SFR is O.lMgy- 1 for the brightest QBCDs, i.e., when 
g ~ -18.5. Figure [H also includes the SFR of BCDs - 
in this case 7 = 1 to acknowledge that the starburst is 
spreadout all over the galaxy. 

Figure [15] shows the difference between nebular and 
stellar oxygen metallicities predicted by the closed-box 
evolution of Class QBCD galaxies. We assume t to be 
the age of the stellar population derived in § [7] Equa- 
tions (USD, ([17]), (Eg]), IJ39J, and ([20]) were used with 
/ = 2 (the square symbols) and / = 0.04 (the times 
symbols). The case / = 2 represents a galaxy- wide 
starburst able to pollute with metals the whole galac- 
tic gas. (/ = 2 assumes the same amount of mass in gas 
as the mass in stars, which is reasonable for low surface 
brightness dwarf galaxies like our QBCDs. According 

6 Note that the SDSS equivalent width corresponds to a spec- 
trum taken at the center of the galaxy ( § [2j ■ Using the integrated 
luminosities to estimate Ha fluxes assumes that the star-forming 
burst observed at the galaxy core extends evenly throughout the 
galaxy. The factor 7 accounts for the case where the starburst 
affects only a faction 7 of the galaxy. 
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Fig. 14.— Star Formation Rate (SFR) vs SDSS g absolute 
magnitude. Both QBCD galaxies (Class 0; the square symbols) 
and BCD galaxies (full set; the times symbols) are included. The 
dashed line represents the average QBCD SFR, i.e., SFR ~ 0.25 • 
(M*/10 10 Mq). The solid line corresponds to a SFR 50 times 
larger. The stellar mass of the QBCD galaxies derived from the 
g magnitude is also included in the scale on top of the figure. 
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Fig. 15. — Difference between nebular and stellar metallicities 
(Z — Z s ) predicted if the chemical evolution of Class QBCDs 
follows a closed-box model. The difference is plotted vs the ab- 
solute g-color magnitude (bottom axis) and vs the stellar mass of 
the galaxy (Af*, top axis). The horizontal solid line corresponds 
to the observed difference, which can be reproduced only if a small 
fraction of the galactic gas is polluted with metals (see the main 
text) . 

to lStavelev-Smith et alj (| 19921) , they have one Mq of Hi 
gas per solar luminosity, which corresponds to / ~ 3.) 
In this case the predicted difference of metallicity is too 
low to account for the observed difference (§0, which is 
represented in Figure [15] as a horizontal solid line. The 
amount of metals produced at the current SFR during 
the age of the starburst is insufficient to effectively con- 
taminate the whole ISM of the galaxies. If, on the other 
hand, the same starburst pollutes only a small fraction 
of the galactic gas (/ = 0.04, or a factor 50 smaller than 
the previous case), then the predicted and the observed 
metallicities agree (Fig. [15j the times symbols). Con- 
sequently, the observed Z — Z s can be explained if the 
gas from which we infer the metallicity represents only a 
small fraction of the total galactic gas. 



Agreement between observed and model metallicities 
can be also reached if the QBCD galaxies had an average 
SFR during t much larger than the present one. The pre- 
dicted closed-box evolution depends on the ratio SFR// 
rather than on SFR and / separately, and a decrease of 
/ is equivalent to an increase of SFR; see equations (|26j) 
and (|28p. However, the required increase of SFR is too 
high to be sensible. It would have to be 50 times larger 
than the observed ones. This level of continuous star for- 
mation activity during the last two Gyr is unreasonably 
high; it is shown as a solid line in Fig. [14] and it corre- 
sponds to <~ 5M Q y _1 for the brightest QBCD galaxies. 
It is larger that the already large SFR observed in our 
BCD galaxies (see Fig. I14| the times symbols). 

So far our argumentation has considered closed-box 
chemical evolution. If the box is opened both infall of 
(low metallicity) gas from the intergalactic medium, and 
outflows of metal rich SNa eject a are to be expected (e.g., 
iGarnettl 120021 : lDalcanto"nl I2007L and references therein) , 
and both processes reduce the metallicity of the gas with 
respect to the predictions of the closed-box model. Then 
it becomes even more difficult explaining the observation 
as a global chemical enrichment. An approximate way 
of considering this open-box evolution using closed-box 
equations consists of using effective yields inferred from 
observations rather t han the yield pre dicted by stellar 
evolution models (e.g. iDalca nton 20075) ■ These effective 
yields are smaller than true yields, and the difference 
increases with decreasing galaxy mass (probably due to 
the decrease of the gravitational potential of the galaxy). 
Actually, the deviations are p articularly large for dwarf 
galax ies, like our QBCDs (e.g., lGarnettj|2002t iDalcantonl 
l2007f ). If the yield y is reduced, then Z — Z s decreases 
too - see equation (|27|) and keep in mind that /i is fixed 
by equation (|2"(5|) . In short, the assumption of closed- 
box chemical evolution do not invalidate our conclusion, 
namely, that the metallicity we infer from the emission 
lines traces only a small fraction of the total galactic gas. 

So far we have put aside the behavior of the BCDs. 
According to Fig. [12] they have the same nebular and 
stellar metallicities within error bars. In addition, the 
stellar populations have ages slightly smaller than 1 Gyr 
(Fig. [13]). These two properties can be easily explained 
if the BCD galaxies are QBCDs experiencing a major 
but short starburst involving fresh well mixed galactic 
gas. In this case, the observed BCD SFR is significantly 
larger than the average SFR during the age of the stellar 
population. Then the chemical evolution model to be 
applied has M ~ M*, and a SFR of the order of the 
SFR of QBCDs, i.e., a model similar to the case labelled 
as 2 M* in Fig. Qjj] The model prediction is Z — Z s at 
least one order of magnitude smaller than the difference 
observed in QBCDs and, therefore, in agreement with the 
lack of metal enrichment observed in BCDs (the stars in 
Fig. H2J. 

9. CONCLUSIONS 

We analyze the metallicity of QBCDs, i.e., galaxies 
that may be Blue Compact Dwarfs (BCDs) during the 
periods of quiescence, when the major starburst charac- 
teristic of BCDs is not so dominant. The QBCD candi- 
dates were selected in Paper I from SDSS/DR6, where we 
also separate a reference sample of BCDs. The metallic- 
ity inferred from emission lines of these QBCDs turned 
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out to exceed the metallicity of the BCDs, an uneasy 
result if BCDs have to descend from QBCDs. Here we 
study whether the metallicity inferred from the emission 
lines of QBCDs may not be representative of the full 
galactic gas, but reveals a local enrichment by recent 
starbursts. In this case the metallicities for the gas and 
for the stars must differ significantly. 

The work is based on SDSS/DR6 spectra, whose 
signal-to-noise ratio is not sufficient to measuring stellar 
metallicities from absorption lines of individual spectra. 
We improve the original signal-to-noise ratio by stacking 
observed spectra that are alike. The grouping of similar 
spectra was carried out by classifying the 21493 QBCD 
galaxies using an automatic k-means classification pro- 
cedure (§ EJ. The algorithm renders a small number 
of types of spectra or classes, with the first ten classes 
containing 90% of all spectra, and the most numerous 
Class having more than 36% of them. As a by-product, 
this classification scheme provides a selective technique 
to identify galaxies in various states within the color se- 
quence. In particular, one of our classes seems to con- 
tain only transition galaxies in the so-called green valley 
(Class 1, § 13. ip . whereas another class includes most of 
the red sequence galaxies (Class 3; § I3.ip . The typical 
Class QBCD galaxies belong to the blue sequence. So 
do BCD galaxies. 

The stellar metallicities have been derived from the 
absorption lines using an ad-hoc procedure which fits 
the average profile of each class with single-stellar pop- 
ulations synthetic spectra based on the stellar library 
MILES (§ S]). We develop our own simple but robust 
tool to get a intuitive control of the errors. Emission 
lines are masked out. The galactic continuum is also 
subtracted for fitting, so that only absorption line fea- 
tures contribute to the measurement. As inferred from 
our Monte-Carlo estimate of the random error budget, a 
direct fit of the full spectrum is good enough to assign 
ages, but it does not provide enough finesse to properly 
distinguish metallicities. Only after overweighting par- 
ticular bandpasses of the stacked spectra (corresponding 
to so me of the Lick indexes defined by iWorthev et al.l 
fl99l . we bring the formal error bars down to reasonable 
limits, below 0.1 dex. Gas metallicities arc obtained from 
emission lines with errors smaller than 0.2 dex (§ [5] and 
App. [A"]) . When the gas and the stellar QBCD metallic- 
ities are compared, gas metallicities turn out to be sys- 
tematically larger than the stellar metallicities by some 
~ 0.35 dex (§[6|). Despite the fact that this difference 
is not far from the formal error bars (actually, it is be- 
low the formal 2 a level; § [6]), we regard it as signifi- 
cant for a number of reasons. First, it is systematic, so 
that the main QBCD classes show it. Second, it is not 
present in BCDs, where stars and gas show the same 
metallicity within error bars (Fig. [TTj) . Third, the ex- 
cess of gas metallicity with respect to stellar metallicity 
is implicit in the luminosity-metalli city relationsh i ps for 
star-form ing galaxies in f erred by iTremonti et alJ (2004, 
gas) and iGallazzi et al.1 (|2005l stars) (see § [6j and the 
dashed line in Fig. [T2"|) . Despite the existence of all these 
supportive arguments, a caveat is in order. The stellar 
metallicity error bars only describe statistical errors, al- 
though systematic errors may dominate the error budget. 
Even if these systematic errors exist and are important, 
they should not modify the conclusions as they would af- 



fect both QBCDs and BCDs in the same way. However, 
one can never discard a source of (unknown and unsus- 
pected) systematic errors affecting QBCDs and BCDs 
differently, which would force us to reconsider the metal- 
licity discrepancies. 

The fraction of QBCD galactic light produced by stars 
augment with the metallicity, so that the fainter the 
emission lines the more metal rich the gas. This result 
reinforces the conjecture that the emission lines come 
from self enriched ISM. The luminosity weighted ages of 
QBCDs span the full range from 1 to 10 Gyr (Fig. [T3]) . 
The most common Class is in the young part of such 
range, with an age below 2 Gyr. The fact that the age 
of the stellar population of BCDs is shorter but not very 
different from the age of QBCDs adds on to this picture. 
If BCDs and QBCDs are basically the same galaxies, but 
the BCDs happen to be in a phase of enhanced star for- 
mation activity, then the underlaying stellar populations 
must have similar properties. Their stellar metallicities 
are similar (§ |6|). The BCD ages are smaller than the 
QBCDs ages, but this can be easily due to the fact that 
their luminosity weighted average ages are reduced by 
the current star burst. 

In principle, the excess of metals in the ionized gas of 
QBCDs, as revealed by their emission lines, may reflect 
the natural enrichment of the ISM produced by succes- 
sive SN ejecta. Emission lines trace the present ISM, 
whereas stars sample it in the past when the metallic- 
ity was lower. The relative enrichment depends on the 
age of the stars, and also on the star formation rate pro- 
viding the SNe. In the case of our QBCDs, these two 
quantities are tightly constrained. We have estimated 
the (mean) age of the starburst, and the (current) star 
formation r ate (SFR) as infe rred from the observed Ha 
luminosity (|Kennicuttl fl998). Using simple closed-box 
chemical evolution models, we argue that given the age 
and the star-formation rate, the observed starburst is not 
sufficient to enrich the full galactic ISM to the observed 
levels. However, age and SFR can be accommodated if 
the enriched galactic gas represents only a small fraction 
of the total gas (~l/50; §[8]). The assumption of closed- 
box evolution does not invalidate the conclusion. 

As we point out in Paper I, QBCD are quiet com- 
mon, representing one out each three dwarf galaxies in 
the local universe. Since they are so common, it is con- 
ceivable that some of their properties are not exclusive of 
the QBCD class, but a global property of dwarf galaxies 
with emission lines. In particular, the bias of metallicity 
inferred from emission lines may be present in all star- 
forming dwarf galaxies, rather than being a feature of our 
particular subset. We plan to explore this potential bias 
using the techniques developed in the paper, namely, the 
comparison between the metallicity estimates based on 
emission lines and absorption lines. Moreover, we plan 
of applying the classification tool in § [3] to find out and 
characterize galaxy spectra corresponding to the vari- 
ous parts of the color sequence. The short green-valley 
phase is particularly interes ting (e.g., iDellenbusch et all 
l2007t iSilverman et alJl2008h . and we came across a sim- 
ple method of identification. 
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APPENDIX 

RANGE OF PHYSICAL CONDITIONS WHERE N2 HOLDS 

The strong-line index abundance estimates are c alibrated against T e methods using particular sets of targets (e.g., 
iDi'az fc Perez-Montero|[2000t iPettini fc Pagelll2004ft . There is no guarantee that such calibrations hold when applied 
to galaxies whose physical conditions deviate from those used for calibration. In particular, if the physical conditions 
prevailing in the H II nebulae of QBCDs and BC Ds are not included among the calibration targets, then the metallicity 
difference between BCDs and QBCDs found by San chez Almeida et ail (|2008t ) may be an artifact of using the strong- 
line N2 method. The physical conditions in the BCD H n regions (newborn) and in the QBCDs H n regions (aging) 
are expected to differ systematically which, coupled with an unfit N2 calibration, may give rise to a false difference. 
Fortunately, such potential bias seems to be harmless in our particular case. The N2 method holds for a range of 
physical conditions broad enough to encompass both BCDs and QBCDs. Figure [TBk shows the residual between the 
abundance determined with the N2 method, [12 + log(0/H)]N2, and the T e method, [12 + log(0/H)]T e , 

A[12 + log(0/H)] = [12 + log(0/H)] N2 - [12 + log(0/H)] Te , (Al) 

for the set of ~310 metal-poor line-emission SDSS/DR3 galaxies analyzed by IIzotov~et al. (200§). The T e based abun- 
dances are directly computed by the authors from state-of-the-art photoio nization modeling, wh ereas the N2 abun- 
dances have been derived from the published line fluxes using the recipe in IPettini fc Pagel (|2004j) . 
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Fig. 16. — (a) Differ ence between the ox ygen abundances d etermined with the N2 and the T e methods, A[12 + log(0/H)], for the set 
of galaxies analyzed by Izotov ct al. (2006, the asterisks) and Shi ct al. (2005, the squares). They arc represented vs the [OIII] electron 
temperature T e . There is no systematic trend for T e < 17000 K. (b) A[12+log(0/H)] vs electron density rt e for the set of galaxies analyzed by 
Shi et al. (2005). No systematic trend, (c) Tight correlation between oxygen abundance and electron temperature for 12 + log(0/H) > 7.5. 
(d) A[12 + log(0/H)] versus 12 + log(0/H). As indicated by the insets, the hashed regions correspond to the range of abundances of BCDs 
(the -45° slanted lines), and the QBCDs (the +45° slanted lines). 

The difference is represented as a function of the [OIII] electron temperature. There is no obvious systematic 
difference between the two methods when 5000 K < T e <17000 K. The hint at systematic deviation when T e > 17000 K 
does n ot seem to affect our estimates (see below) . Figure [T5b also includes the difference computed by I Shi et all 
(|2005f ) 7 for a set of ~ 70 BCD galaxies (the squares), which we include because these authors also provide the second 
independen t thermodynamic parameter, namely, the electron density n e . Figure [16b shows the variation of the residual 
with n e for I Shi et all (|2005l ) galaxies, which evidences no systematic trend even though the large range of electron 
densities that are involved; 10 cm" 3 < n e < 10 4 cm" 3 . Finally, Fig. fTBH shows the residuals vs the metallicity. 
Again residuals are independent of metallicity except at the low end when 12 + log(0/H) < 7.5. The strong-line 
methods work because the see mingly indepen dent variables (T e and n e ) are actually related to 12 + log(0/H) as a 
single parameter family (e.g., iStasih ska 2004 § 2.1.2). This is the case for the dependence of T e on 12 + log(0/H); 
see Fig. IT6"b . This tight correlation breaks down at 12 + log(0/H) < 7.5. Due to the correlation, the low abundance 
regime coincides with the high temperature points which seems to deviate from the nominal law as hinted in Fig. 116b . 
This potential problem at low abundances and high temperatures should not affect our sample of BCDs and QBCDs, 
whose abundances correspond to the well behaved part of the N2 calibration: see Fig. [16H . where the range of BCD 
and QBCD abundances are represented as hashed regions. Moreover, according to Figs. [16k and [T6b. the range of 
T e and n e where the N2 calibrations ho ld covers the full the range of thermodynamic parameters to be expected in 
H II regions (e.g.. lDvson fc Francol l2000). These two facts indicate that the N2 abundance estimate is valid for both 
BCDs and QBCDs, even though their H II regions may have different physical conditions. Then the relative difference 
between the abundances of BCDs and QBCDs are not an artifact o f using strong-line empirical methods. The rms 
variations of A[12 + log(Q/ H)1 are 0.16 dex and 0.13 dex, for lShi et al.l (|20 05l) galaxies and llzotov et al.l (|2006h galaxies, 
respectively. In the case of llzotov et al.l (|2006l ). we only consider 12 + log(0/H) > 7.5. 



7 In this case the authors compute both the N2 and the T e abundance estimates. They differ systematically by 0.2 dex. a constant term 
that we have removed when plotting since it does not affect relative differences between N2 abundances. 



